Quantum emulsion: a glassy phase of bosonic mixtures in optical lattices 
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We numerically investigate mixtures of two interacting bosonic species with unequal parameters 
in one-dimensional optical lattices. In large parameter regions full phase segregation is seen to 
minimize the energy of the system, but the true ground state is masked by an exponentially large 
number of metastable states characterized by microscopic phase separation. The ensemble of these 
quantum emulsion states, reminiscent of emulsions of immiscible fluids, has macroscopic properties 
analogous to those of a Bose glass, namely a finite compressibility in absence of superfluidity. Their 
metastability is probed by extensive quantum Monte Carlo simulations generating a rich correlated 
stochastic dynamics. The tuning of the repulsion of one of the two species via a Feshbach resonance 
drives the system through a quantum phase transition to the superfluid state. 
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Trapped ultracold atoms in optical lattices (OLs) of- 
fer unprecendented opportunities of studying the quan- 
tum behavior of correlated matter at a very fundamental 
level, namely by experimentally implementing Hamilto- 
nians which would have been regarded in the past as toy 
models [l|. The introduction of quenchend randomness 
widely enlarges the range of quantum phases which can 
be explored in such systems. In the weakly interacting 
regime the presence of disorder promises to realize Ander- 
son localization of coherent matter waves [2] ; in presence 
of interaction between bosons, condensate fragmentation 
into macroscopically many localized states or Anderson 
localization of collective modes leads to the insulating 
Bose-glass phase Yet a challenging issue is how to 
introduce model disorder in the system. In general all 
disorder potentials realizable experimentally show some 
form of spatial correlation: if the disorder correlation 
length is much larger than the typical correlation length 
of the particles, the disorder imposes a simple density 
modulation fully captured by a local density approxima- 
tion 0] , implying that no new phase is established in the 
system. This is the conclusion of the seminal experiments 
carried out with speckle-laser potentials 0, [B| . An alter- 
native way of producing disorder optically is by loading 
the atoms in optical quasicrystals, for which an insulat- 
ing behavior with unconventional spectral properties has 
been reported recently [6[. 

In this paper we explore an alternative proposal to in- 
troduce disorder in a system of cold bosons in an OL. A 
second species of b bosons is added to the system, having 
a much smaller hopping amplitude than the primary a 
bosons and therefore producing an effective quasi-static 
disorder potential if prepared in random density configu- 
rations or in a quantum superposition thereof 0, Q ■ The 
original proposal of Ref. 8 involved the separate prepa- 
ration of the two bosonic species followed by the sudden 
onset of the interspecies interaction, possibly leading to 
out-of-equilibrium localization of the fast-moving bosons. 
Here we explore the simpler case in which the two species 
are adiabatically prepared together, aiming at their joint 



ground state at T = 0. We base our analysis on quan- 
tum Monte Carlo simulations of a one-dimensional sys- 
tem, targeting the ground state either through thermal 
annealing, i.e. by slowly cooling the system to T = 0, and 
through quantum annealing, i.e. by a slow change of the 
Hamiltonian parameters at low temperature. Our cen- 
tral finding is that, for weak interactions between the a 
bosons, the system displays an exponentially big number 
of metastable states which make the true ground state 
of the system essentially unreachable. While the ground 
state would show perfect phase segregation between the 
two species, in the metastable states the two species are 
fragmented into small droplets with short-range density- 
density correlations, implying phase separation over the 
length scale of a few lattice spacings (depending on the 
annealing rate) . These quantum emulsion states have an 
energy that depends roughly linearly on the phase inter- 
face, reminiscent of the behavior of metastable emulsions 
of immiscible fluids 9] . Both upon thermal or quantum 
annealing, a significant deviation from the true ground 
state energy is observed (residual energy), which decays 
as a small power of the annealing rate; this behavior is 
analogous to that of spin glasses [lj| [Hj . The density 
arrangements of the b bosons in the metastable ensemble 
realize the statistics of a weakly correlated disordered 
potential. Consistently, the global properties of the a 
bosons in the quantum emulsion states are analogous to 
those of a Bose glass, namely gapless excitations and a 
finite compressibility in absence of superfluidity 0] . 

We describe the bosonic mixture in a one-dimensional 
OL via a two-flavor Bose-Hubbard model [l2| 

tt = J2[~ Ja ( a * a l+i + h - c ) ~ Jb ( b i b l+i + h -' 

a 

iV*( n <M _ 1 ) + ~jT n bA n b,i - 1) + U ab n a jnb,i .(1) 

where J a n,) is the hopping amplitude, U a (h) is the on- 
site intra-species repulsion, and U a b is the inter-species 
repulsion. Experimentally, a bosonic mixture of, e.g., 
87 Wo or of Na in two different hyperfine states [IH loaded 
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FIG. 1: (Color online) (a) Correlation between the energy 
and the phase interface in the metastable ensemble. The 
model parameters are U a = J a , Ub = Uab = 5J a , Jt = Ja/5, 
(n a ) = 1 and rib = 1/2. Each point is obtained after 1000 
thermalization MCS and 2000 measurement MCS after cool- 
ing at a rate r — 2 * 10~ 4 . (b) Density-density correlation 
function of the b bosons, Cb(r) = {(rib,i — fib)(nb,i+ r — fib)). 
(c) Two density snapshots of a long QMC simulation. 



in spatially anisotropic OLs realizes the above Hamilto- 
nian (neglecting the trapping potential). Due to the dif- 
ferent dipolar coupling of the two species to the OL, the 
hopping amplitudes J a and Jj, can be significantly differ- 
ent from each other. In what follows we take the case 
of slow b bosons, Jh = J a /5- The depth of the OL fixes 
the ratios U a b/J a = Ub/J a = 5 [DJ. Finally, the appli- 
cation of a magnetic field in the proximity of a Fcshbach 
resonance allows to tune U a . This enables to explore the 
various regimes of strong vs. weak interaction for the 
mobile bosons, which in presence of quenched disorder 
correspond to different localized phases competing with 
superfluidity [IBj . 

We study the Hamiltonian Eq. ([1]) making use of the 
Stochastic Series Expansion (SSE) quantum Monte Carlo 
based on the directed- loop algorithm [l6| with optimized 
transition probabilities [17| . To efficiently reproduce the 
correlations between the two bosonic species we also in- 
troduce a novel double-worm algorithm which allows for 
a simultaneous correlated update of worldlines of both 
species [l8|. One Monte Carlo step (MCS) is composed 
of a diagonal update and as many single- and double- 
loop updates as to visit on average all SSE vertices twice. 
Our simulations are performed in a mixed ensemble, 
namely we fix the density fib = 1/2 of b bosons 14j |. 
while a chemical potential term — fi a n a ,i is added to 



adjust the average filling (n a ) to a fixed value. Through- 
out the paper we choose \x a such that (n a ) ~ 1 (within 
less that 2%). Chains with periodic boundary conditions 
and sizes up to L = 96 are investigated, at inverse tem- 
peratures (3, J a = L/2 or higher (no significant difference 
is appreciated with [3J a — L). We allow for up to 5 (3) 
a (b) bosons per site to reproduce the softcore limit. 

A first empirical observation concerns the difficulty of 
QMC in equilibrating the system in the case U a < U a b, 



despite the non-local algorithms used for the QMC up- 
date. This suggests immediately that many metastable 
states are present in the system, and even a non-local 
stochastic dynamics, as the one generated in the QMC 
update supplied with multi-cluster moves, is unable to 
efficiently escape from local energy minima. At the 
same time, the metastable states are very robust, resist- 
ing essentially unchanged to more than 10 5 MCS (see 
Fig. [TJ). This is strongly reminiscent of glassy sys- 
tems as, e.g., spin glasses, for which the phase space is 
essentially fragmented into metastable "valleys" which 
are not efficiently connected to each other by the dy- 
namics of the system. A fundamental ingredient which 
defines the class of observed metastable states is the 
annealing procedure followed while preparing the sys- 
tem. Given that SSE is intrinsically a finite-temperature 
method, we adopt a thermal annealing scheme by start- 
ing at a temperature T max = 2 J a [l9[ and linearly cool- 
ing the system in steps of AT = 0.02J a down to the 
physical T = 0. At each temperature step we wait M 
MCS, which define the cooling rate r as r = AT/M. 
Fig. [TJ shows the energy E = (Ti) of the metastable 
states, reached by repeatedly cooling down to T = AT 
at a rate r = 2 x 10 -4 , as a function of the phase 
interface A = £\ (9(n a ,i)0(n^ i+1 ) + 0(n bti )9(n a .i+i)) 
(#=Heaviside function), giving the average number of 
nearest neighboring site pairs having different species on 
the two sites, and essentially counting the number of 
droplets. There is a clear correlation between energy and 
phase interface, indicating that the energy is minimized 
for minimal A, namely for phase segregation [20L 21 1. 



Yet for all the metastable states the phase interface is 
a significant fraction of the system volume, which means 
that these states are microscopic emulsions of the two 
bosonic species, with short-ranged density-density corre- 
lations (Fig. Q]). 

It is interesting to observe that the rescaled data for 
different system sizes all collapse onto the same, ap- 
proximately linear dependence of E on A, whose slope 
gives the surface tension of the microscopic emulsion. 
This surface tension has an exquisite quantum mechan- 
ical origin (hence the name of quantum emulsion). In 
fact, while the surface tension in classical emulsions 
comes from long-range attractive interactions between 
the atoms/molecules [9j, in our case the interactions are 
all repulsive and on-site only. In absence of the quantum 
kinetic term, the energy does not depend at all on the 
phase interface, but only on the phase overlap. Therefore 
the surface tension stems from the attempt to minimize 
the quantum zero-point kinetic energy. 

The almost univocal dependence of E on A indi- 
cates that spatial permutations of the emulsion droplets 
leaving the phase interface unchanged produce quasi- 
degenerate states. The collapse of the data for differ- 
ent sizes shows moreover that the typical droplet size is 
only weakly dependent on the system size. Therefore the 
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number of droplets in the emulsion scales linearly with 
the system size, and so the number of droplet permuta- 
tions scales exponentially. This implies that the number 
of quasi-degenerate emulsion states scales also exponen- 
tially. Hence the strong tendency of the system to get 
trapped in the metastable states. 
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FIG. 2: (Color online) (a-b) Energy vs. cooling rate for 
thermal annealing (Hamiltonian parameters as in Fig. [Tj. 
AT = 0.05J a for L = 48 and 0.02J a for L = 96. (c) Time 
schedule in the quantum annealing procedure; in the sim- 
ulation the time variable is changed in steps of 5 * 10~ 2 . 
MF=magnetic field, (d) Energy vs. annealing rate for quan- 
tum annealing. 

Fig. [5] shows the dependence of the average energy of 
the metastable ensemble at T = AT (see Fig. caption) as 
a function of the cooling rate r. Averaging over repeated 
cooling cycles is made, with ~ 200 takes for the fastest 
cooling rates and ~ 20 takes for the slowest ones. It is 
clear that the final energy systematically deviates from 
the true equilibrium value and it depends on r, although 
the dependence is extremely slow. A fit to a power-law 
dependence of the kind E(r)/L = eo + b gives cj) = 
0.18(3) for L = 48 and = 0.09(5) for L = 96; this might 
suggest that <fi — > for L — > oo, leaving out a logarithmic 
dependence only which is generally expected in presence 
of low-energy metastable states [111, |22j . A similar slow 
dependence on r is also observed for the average phase 
interface, as expected from the linear relation with the 
energy. Therefore the kinds of quantum emulsion reached 
at different cooling rates are weakly dependent on r itself, 
and the behavior described before for r = 2 * 10~ 4 is 
generic for a broad range of r values. 

Although thermal annealing is quite convenient for 
QMC simulations, it is not necessarily the most relevant 
annealing procedure for experiments. In fact, a typical 
OL experiment is in principle performed at fixed, very 
low temperature, changing the Hamiltonian of the bosons 
from that of a trapped weakly interacting gas to that of a 



strongly interacting one in the OL. Therefore the picture 
of so-called quantum annealing 11, [23[ is more appro- 



priate for the experimental preparation of the system. 
We simulate an incoherent two-step quantum annealing 
process in which first the ratios J a /U a (h) and Jb/J a are 
linearly decreased, mantaining U a = Ub = U a b (OL ramp- 
up) , and then U a is decreased by tuning an applied mag- 
netic field close to a Feshbach resonance. The whole an- 
nealing is performed with M MCS, defining the annealing 
rate tq = 1/M. We observe that also upon quantum an- 
nealing the system remains far from the true equilibrium; 
we cannot fit E vs. tq neither with a power law (as in 
the thermal case) nor with a logarithm, which suggests 
that even for the slowest annealing rates we are far from 
the asymptotic behavior for rQ — > 0. 

Having shown the robust features of the ensemble of 
metastable emulsion states, we argue that this ensem- 
ble is experimentally relevant. The existence of expo- 
nentially many metastable states implies the presence of 
many avoided level crossings between such states upon 
an arbitrary evolution of the Hamiltonian parameters 
toward their final value. The gaps at the avoided level 
crossings are exponentially small in the ratios J a (b) /U a (b) > 
given that the metastable states typically differ from each 
other by a macroscopic rearrangement of the density dis- 
tribution. The evolution of the system's state during the 
experimental transformation of the Hamiltonian param- 
eters will then exhibit a cascade of Landau-Zener (LZ) 
processes, so that the final state |\E f ) of the system will 
be a large quantum superposition \^>) = X)v c V'l^) °f 
different metastable states \tp) |2~ij ]. 




FIG. 3: (Color online) Phases of the a bosons upon tuning 
the repulsion U a - Upper panel: superfluid fraction p Sja and 
compressibility « a . Lower panel: coherent fraction N a (k — 0) 
and inverse correlation length £J . In the inset: sketch of the 
interaction-assisted resonant tunneling. 

The above argument motivates us to regard the 
metastable emulsion states as the experimentally relevant 
ensemble for a given annealing procedure: incoherent 
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averaging over this ensemble essentially corresponds to 
measurements in the experimental final state of the sys- 
tem 24|. The averaging over the emulsion states reveals 



an out-of-equilibrium quantum emulsion phase, with ex- 
otic macroscopic properties. Fig.[3]shows the global prop- 
erties of the a bosons as a function of the tunable U a /J a 
repulsion upon averaging over 100-200 thermal anneal- 
ing cycles at a rate r = 2 * 10~ 4 . It is observed that 
for U a /J a < 3.5 the a bosons are in a quantum emul- 
sion phase with zero superfluid fraction Ps.a (estimated 
through winding number fluctuations [161 l25j|) but finite 
global compressibility n a = (3([J2i( n a,i- (n>a))] 2 ) / L- The 
absence of superfluidity is clearly due to the fragmenta- 
tion of the a bosons into incoherent metastable droplets. 
Such fragmentation can be regarded effectively as a phe- 
nomenon of localization due to the strong repulsion with 
the b bosons, which are also fragmented into weakly cor- 
related droplets. In this respect, the quantum emulsion 
phase behaves as a metastable Bose glass: although not 
superfluid, this phase is compressible because arbitrar- 
ily large droplets of a bosons can appear which admit 
gapless particle number fluctuations. 

Upon increasing the on-site repulsion U a /J a , the lo- 
calization of a bosons can be overcome by a phenomenon 
of interaction-assisted resonant tunneling. An effective 
hopping w Ja/{U a b — U a ) (to 2 nd order perturbation 
theory) brings an a boson from a doubly occupied site 
to a singly occupied one across an intermediate site oc- 
cupied by a 6 boson (see Fig. [3]). Close to the reso- 
nance condition U a w U a b this process becomes very 
effective, and the coherent tunneling of a bosons from 
one droplet to the neighboring one leads to the onset of 
superfluidity. Roughly speaking, when ijj C 1 a frac- 
tion 1-nb of a bosons screens the potential created by 
the 6-boson droplets, and it assists the resonant tun- 
neling of the remaining fraction rib, which then leads to 
Ps,a ~ fib- The phase of the system with superfluidity of 
the a bosons still shows a significant metastable behav- 
ior: the spatial arrangement of the 6-boson droplets is 
in fact quasi-static in the QMC simulation, namely the 
simulation is strongly non-ergodic in that translational 
invariance remains broken over more than 10 5 MCS. It is 
suggestive to picture the b bosons and the normal frac- 
tion of the a bosons as frozen in a metastable quantum 
emulsion, whereas the superfluid fraction of a bosons is 
able to coherently tunnel through the emulsion and to 
ergodically visit the entire lattice. Concerning experi- 
mentally relevant quantities, the onset of superfluidity 
is marked by a dramatic enhancement of the coherence 
peak N a (k = 0) = J2ij( a l a j) I 'L, and by a strong de- 
crease of its width (quantified by the inverse phase cor- 
relation length through second- moment estimation). 
Upon further increasing U a /J a the unit-filled a bosons 
are eventually driven to a Mott insulating state, forming 
an almost uniform background for the b bosons which 
delocalize into a superfluid state. 



In conclusion, we have reported a strongly metastable 
behavior of repulsive bosonic mixtures in one- 
dimensional optical lattices. An out-of-equilibrium 
quantum emulsion phase emerges from QMC simula- 
tions, and it is found to be very robust to thermal as well 
as quantum annealing. This phase has a pronounced 
glassy behavior, with a significant residual energy 
slowly decaying with the annealing rate. All these 
ingredients suggest that this phase be extremely relevant 
for ongoing experiments on bosonic mixtures as well as 
for one-dimensional fermion-boson mixtures in optical 
lattices. 
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